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1. Introduction 

The QCD phase diagram contains several kinds of states of matter on the pL-T plane, where T 
is temperature and pL is quark chemical potential. So far, attentions have mostly been paid to a high 
temperature and low density region, stimulated by heavy ion experiments. It is also interesting to 
study the low temperature region of the QCD phase diagram, such as nuclear matter, high density 
states inside neutrons stars, and color superconducting matter. 

Lattice simulations suffer from infamous sign problem at nonzero quark chemical potentials. 
In addition, other difficulties arise at low temperatures. Numerical costs increase with A',, and 
the fermion determinant becomes less sensitive to ji. at low T and small jX. The sign problem 
becomes severe at jLt ~ mji/2, and Monte Carlo simulations are difficult for the region. This causes 
a problem called early onset problem of the quark number density between mji/2 < pL < Mn/3. 
For instance, it is unclear how nuclei or nuclear matter are spontaneously formed from quarks at 
low temperatures and low densities. 

Recently, we have addressed the properties of QCD at low temperatures [|l]], using a reduction 
formula of the fermion determinant ^, 0, ^, ^. The reduction formula leads to a reduced 
matrix, which is a product of block-matrices describing hops of quarks from a time slice to the 
next time slice. We found a Nt scaling law of the eigenvalues of the reduced matrix, where Nt is a 
temporal lattice size. The Nt scaling law was used to take the low temperature limit of the fermion 
determinant. In |[l]], we obtained a limit of the fermion determinant at low T and small jj., which is 
independent of jj.. The result is consistent with studies in ||], and provides the information for 
the properties of QCD at low T and small jj, based on the lattice QCD simulation. In addition, we 
obtained a low T and large /i limit of the fermion determinant. Such a study would provide a clue 
to understanding the finite density QCD at low temperatures. 

In this paper, we will review the idea to take the low temperature limit. This work is a series of 
our studies of finite density lattice QCD with a clover-improved Wilson fermion [^, ^, |l0[ H), 



2. Reduction Formula 



We begin with the clover-improved Wilson fermion defined by 



i=l 



K 



e+^''{r-Y4)U4{x)d^,^^+i + e'^''{r+Y4)Ul{x')d^,^^_^ -KCswd.y £ O^vF^y, (2.1) 



where K and r are the hopping parameter and Wilson parameter, respectively. We denote diagonal 
elements of A with regard to time as B, 

3 

B = - K-£ (r - Y)Uiix)5,.,^j+ ir+y^ulix')5^.,_. - kCsw5,,,' £ o^yF^y, (2.2a) 

1=1 ii<v 

and introduce two types of block-matrices 

ai = B"'''^''{x,y,ti) r^" - 2k r^''5"''5{x-y), (2.2b) 
A- = B^'''^^ix,y,ti) r^''Uf{y,ti) - 2k r^_' 5{x-y)Ut{y,ti). (2.2c) 
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r± = (rib 74) /2 are projection operators in case of r = 1. Using the reduction formula, the fermion 
determinant is given by 



detA(At) = r'^'"'/'Codet(^ + G), 
where t, = exp(-/x/r), A^-ed = ANcN^,, and 

!2 = («f'/3i)---(a^//3,v,), 
Co = ^ndet(a;)^ . 

Q and Co are independent of jU. Denoting the eigenvalues of 2 by A„, we obtain 

A'red 



detA(At)=Cor^'=''/'n(^" + ^)' 



(2.3a) 

(2.3b) 
(2.3c) 



(2.4) 



which describes the /X -dependence of the fermion determinant. Having the eigenvalues A provides 
values of det A for any . 
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Figure 1: Eigenvalues for Nt = 4- and A', = 8. Top panels : the scatter plot of eigenvalues on the complex 
X plane. Bottom panels : the density of the magnitude of the eigenvalues. K was fixed with a LCP of 



mps/my =0.8 |12|. /3 = 1.86, which corresponds to Tc at = at M = 4. The top panels show eigenvalues 
with |A„| > 1. 

Figure [T] show the Nt dependence of the eigenvalues. Top panels show the scatter plot of 
the eigenvalues A for Nt = 4(left panel) and N = 8(right panel); each panel was obtained from 
a configuration. The two results were obtained with the same values of the parameters except 
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for A^,. The magnitude of the eigenvalues strongly depends on A^,; | A | = 10^ ~ 10'^ for Nt = 4, and 
I A I = 10"* ~ 10^ for Nt = 8. Note that the top panels show eigenvalues with | A„ | > 1 . Bottom panels 
show the density of the magnitude of the eigenvalues; the horizontal axis is In |A | and (In |A D/A'^f in 
the left and right panels, respectively. The eigenvalue density is independent of A', if it is described 
as a function of (ln|A|)/A',. This agreement suggests that there exist A', -independent quantities 
/„ = In |A„|/A'f, and that the eigenvalues can be described as |A„| ~ l^' . This provides a scaling law 
of the eigenvalues with regard to Nt- 

Figure |T] also shows two spectral-properties of the reduced matrix Q. The eigenvalue density 
is symmetric with regard to In | A | =0, which results from the fact that two eigenvalues A and 1/A* 
form a pair. The eigenvalue density has a gap near |A| = 1, where no eigenvalue exist. The gap is 
related to the mass of the pion m^j. Gibbs obtained a relation ^ 



arrijr = max In I At I 



(2.5) 



which is given for single configuration. Fodor, Szabo and Toth obtained a modified expression []13|], 

(2.6) 

which is defined for an ensemble average. Although the two expressions differ, it is expected that 



arriji = lim 



they agree with at low temperatures [13] 



2.1 Low temperature limit of the fermion determinant 

The A'f scaling law of the eigenvalues was obtained from the comparison of the eigenvalue 
distribution for Nt = 4 and N = S. It is desired to confirm the results for larger A';. However, the 
eigenvalue calculation of Q suffers from numerical errors for large A^,, because min(A)/max(A) 
exponentially decreases with A'f. It is noteworthy that the N scaling law is supported by two facts; 



Q is defined as the product of A^, block-matrices by Eq. ( [2.3bD , and 2 is a generalization of the 
Polyakov loop In this paper, we employ the A'^f scaling law to consider the low temperature 
limit A'f — > oo. We consider a fixed lattice spacing a, and assume a fixed lattice volume. 
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Figure 2: A schematic figure for the relation between ^ and A in the low temperature limit in the histogram 
of In \X\/Nt. Dotted lines denote the behavior of — /x/T = ln^ with increasing jj.. 
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The A^, scaling law provides the Nt dependence of detA. We denote an eigenvalue A„, (|A„| > 1) 
as Xn = l^'e'^" > 1), and its counter-pait as 1/A* =/,7^'e'^". detAis given as 

A'red/2 iV,-ed/2 

detA = Co^ -^-'/^ n (^"''"'^' + 'n'e'^" ) fl (^"^"'^' + ^n^'e'^" ) • (2-V) 



The first term in Eq. ( [2.7| ) is, in the limit Nt — > °°, reduced to 

A'red/2 

Yl {e-^"''' + l^'e'^")r^ n Ce'^". (2.8) 
(1=1 )i=i 

Here, the next to leading term is proportional to the spatial volume, which vanishes if we take 
Nt with a fixed N^. The hmit of the second term, (e^^"'^' + l^'^'e'^"), depends on ^u; e^^"^' » 
for small jj., and e^^"^' <ti l^^'e'^" for large }JL. Since the eigenvalues have a finite range of 
the distribution, simple expressions aie obtained at enough small and large quark chemical poten- 
tials. The situation is illustrated in Fig. ^ 

If exp(/ifl') ^ min(/„), detA is reduced to 

A'red/2 

detA = CoY\Xn, ne > 1} (2.9) 

n=l 



which gives the low T and small /i hmit of detA. Equation ( [2.9[ ) is real, because Co and the product 
of large eigenvalues are both real. Furthermore, detA(/i) is independent of jJ., which is consistent 
with the 75 hermiticity of detA at /i = 0. The condition reads /it? = aniji/l owing to Eq. (^). 
Thus, Eq. ( ^ ) leads to the jJ. independence of detA up to /i = mji/2. The same result was already 
obtained in ||7|, |8|]. Our study explains the phenomenon using lattice QCD simulations. 
If exp(;Ufl') » max(/„), detA is reduced to 

detA = (^-^-■'/^CodetQ 

N, 

= ^-^'-^/^Yldet{B"'^^'''{x,y,ti) rl" -2Kr'i''5{x-y)), (2.10) 
which gives the low T ans large /i limit of detA. With Eq. ( |2.10| ), the partition function is given by 
limZcci^iJ) = e^^f^^-^'^/^ J m (detA(At)|r^of^>"^°, (2.11a) 

detA(M)|r^o = fidet {B^'^^^^yJi) r^' -2k r^^^" 5 {x - y)) , (2.11b) 

(=1 

where the low T limit is considered only for the quark, and the gluon part remains unchanged. 
Only the factor e\p{2NfNcN^ ix/T) includes jj., which gives the quai-k number density, 

{n)=2NfNc, (lattice unit) (2.12a) 



which is the maximum number of the quarks inside a finite lattice. This implies that Eq. ( [2.11[ ) is 
obtained if all the states are occupied by quarks. 

The maximum number of the quarks is a consequence of the finite lattice size, and a value of 
/i at which Eq. ( 2.10| ) is valid would increase with N^. Therefore, Eq. (2.11) may differ from a 
physical limit of QCD at low temperatures and high densities. However, it is worth considering, 
because Eq. (2.11 ) is real, and may be used for Monte Carlo sampling at nonzero /i. 
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3. Volume dependence of the eigenvalues 
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Figure 3: The volume dependence of the eigenvalues. Left ; absolute value, right: phase. K was fixed with 
aLCP of OTPS /mv = 0.8. /3 = 1.86. 
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Figure 4: Volume dependence of the eigenvalues. 

In the previous section, we have considered the low temperature limit with a fixed lattice 
volume. In this section, we consider the volume dependence of the eigenvalues. Figure |3| shows 
the density of the eigenvalues for 8^ x 4 and 10^ x 4, and Fig. Q shows the eigen spectrum. The 
eigenvalue density is insensitive to the spatial size Ns both in the magnitude (left panel) and phase 
(right panel). The maximum eigenvalue and the size of the gap are also insensitive to Ns- 

Insensitivity of the eigenvalue density to A'^ provides the Ns dependence of the fermion deter- 
minant. The fermion determinant is given in the following spectral representation, 

detA(Ai) =Coexp (^NcV,ii/T + ANcV^ j t/Ap(A)ln(A + <^)^ , (3.1) 

where p(A) is the eigenvalue density on the complex A-plane. The A-integral, / <iAp(A)ln(A + <^), 
is insensitive to A'^v, because p(A) is insensitive to A^^- This means that detA(/i) is proportional 
to exp(Vi), leading to a well-known fact that the sign problem becomes severe for large lattice 
volume Q. 
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4. Summary 

We have discussed the low temperature Umit of the fermion determinant by using the reduction 
formula of the fermion determinant. We obtained two types of the low temperature limit ; one is for 
small quark chemical potentials, and the other is for large quark chemical potentials. The fermion 
determinant is independent of the quark chemical potential in the low density limit, and holds for 
/I < nij^/l. It has the maximum quark number at the high density limit. Both them are real. In 
particular, the high density limit may be used for generating gauge configurations at nonzero quark 
chemical potentials. 

The results were based on the lattice simulations for 8^ x 4, S'* and 10^ x 4. Further studies 
are important to confirm the A', scaling law and the relation between the eigenvalue gap and pion 
mass, and to consider thermodynamical limit. 

This work was supported by Grants-in-Aid for Scientific Research 20340055, 20105003, 
23654092 and 20. The simulation was performed on NEC SX-8R at RCNP, NEC SX-9 at CMC, 
Osaka University. 
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